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Abstract 

We investigate the liquid-gas phase transition of dense matter in supernova explosion 
by the relativistic mean field approach and fragment based statistical model. The 
boiling temperature is found to be high (T^ou > 0.7 MeV for p B > 10~ 7 fm -3 ), 
and adiabatic paths are shown to go across the boundary of coexisting region even 
with high entropy. This suggests that materials experienced phase transition can be 
ejected to outside. We calculated fragment mass and isotope distribution around the 
boiling point. We found that heavy elements at the iron, the first, second, and third 
peaks of r-process are abundantly formed at p B = 1CP 7 , 1CP 5 , 10 -3 and 10 -2 fm -3 , 
respectively. 
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1 Introduction 



It is generally believed that there exist several phases in nuclear matter. 
Among the phase transitions between these phases, the nuclear liquid-gas 
phase transition has been extensively studied in these three decades [1]. It 
takes place in relatively cold (T boi i = (5-8) MeV) and less dense (p B ~ Po/3) 
nuclear matter, and it causes multifragmentation in heavy-ion collisions. When 
the expanding nuclear matter cools down and goes across the boundary of co- 
existing region, it becomes unstable against small fluctuations of density or 
np asymmetry, then various fragments are abundantly formed almost simul- 
taneously. Especially at around the critical point, fragment distribution is 
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expected to follow the power law [2], Yf oc A~ T , which is one of the charac- 
teristic features of critical phenomena. Recent theoretical model studies [3-7] 
have shown that it is very difficult to describe this fragment distribution in 
a picture of sequential binary decays of one big compound nucleus, which 
has been successfully applied to the decay of nuclei at low excitation. This 
finding suggests that it is necessary to consider statistical ensemble of various 
fragment configurations rather than one dominant configuration in describing 
fragment formation at around the boundary of coexisting region. 

In the universe, the temperature and density of this liquid-gas phase transi- 
tion would be probed during supernova explosion. In the collapse and bounce 
stages of supernova explosion, the density and temperature are high enough 
to keep statistical equilibrium [8,9]. At baryon densities of p B > 1CT 5 fm~ 3 , 
since the density is too high for neutrinos to escape, neutrinos are trapped 
in dense matter. This leads to an approximate conservation of lepton fraction 
Yl = L/B and entropy per baryon S/B, where L and B denote the lepton and 
baryon number, respectively. After the core bounce, supernova matter, which 
is composed of nucleons and leptons, expands and cools down. As the baryon 
density and the temperature decrease, charged particle reactions become in- 
sufficient and the chemical equilibrium ceases to hold, namely the system 
freezes out at this point. If the supernova matter goes across the boundary of 
coexisting region and the boiling point T^i of the liquid-gas phase transition 
is higher than the freeze-out temperature Tf Q , this matter will dissolve into 
fragments and form various nuclei in a critical manner. It further keeps equi- 
librium and expands to the freeze-out point. The statistical distribution of 
fragments at freeze-out would provide the initial condition for following nucle- 
osynthesis such as the r-process. (See following references on r-process [10,11] 
and references therein.) 

The importance of the nuclear liquid-gas phase transition in supernova ex- 
plosion was already noticed and extensively studied before [12]. However, there 
are two more points which we should consider further. First, the main inter- 
est in the previous works was limited to the modification of the equation of 
state (EoS). The nuclear distribution as an initial condition for the r-process 
was not studied well. Secondly, in constructing the EoS of supernova matter, 
the mean field treatment was applied in which one assumed one kind of large 
nucleus surrounded by nucleon and alpha gas [13,14]. At temperature much 
above or below the boiling point, fragment mass distribution is narrow and 
the one species approximation works well. However, since fluctuation domi- 
nates at around the boiling point, it is necessary to take account of fragment 
mass and isotope distribution. This distribution of fragments can modify the 
following r-process nucleosynthesis provided that the freeze-out point is not 
far from the boiling point. 

In this work, we study nuclear fragment formation through the nuclear 
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liquid-gas phase transition during supernova explosion. This process may lead 
to the production of medium mass nuclei as seed elements and serve as a pre- 
process of the r-process. We call this process as LG process [15]. In order to 
pursue this possibility quantitatively, it is necessary to determine the liquid- 
gas coexisting region. We find that the liquid-gas coexisting region extends 
down to very low density keeping the boiling point around T^u ~ 1 MeV in 
a two-phase coexistence treatment of EoS with the Relativistic Mean Field 
(RMF) model [13,16-19]. In supernova explosion, it can happen that material 
with S/B > 10 is ejected to outside [20]. Adiabatic paths of ejecta are found to 
go through the calculated liquid-gas coexisting region even with high entropy. 
Having this finding of the passage through the coexisting region, we investigate 
the fragment distribution at around the boiling point in a statistical models 
of fragments [4-7] , referred to as the Nuclear Statistical Equilibrium (NSE) in 
astrophysics. We show that heavy elements around the first, second, and third 
peaks of r-process are abundantly formed at p B = 10~ 5 , 10~ 3 and 10~ 2 fm~ 3 , 
respectively, with temperatures around and just below Tb u in NSE. Further- 
more, the isotope distribution of these elements are also well described in this 
model. We find that it is important to take account of the Coulomb energy re- 
duction from the screening by electrons in supernova matter. Although nuclei 
formed at high densities p B ~ 10~ 2 fm~ 3 having very small entropy at around 
Tbou is not likely ejected to outside, heavy elements up to the r-process third 
peak are already formed statistically at these densities. 

This paper is organized as follows. We describe the treatment of two-phase 
coexistence with RMF model in Sec. 2. The liquid-gas coexistence is shown in 
the (p B ,T) diagram. We demonstrate that it would be possible for a part of 
ejecta in supernova explosion to experience the liquid-gas coexisting region. 
The effects of liquid-gas coexistence on EoS, proton fraction, and adiabatic 
path are also studied. In Sec. 3, we describe the nuclear statistical model 
of fragments at equilibrium (NSE) to study the production of elements. We 
take into account the Coulomb energy modification from electron screening. 
We evaluate the fragment distribution at around the boiling point and in the 
coexisting region within this statistical model. We found that fragments are 
formed abundantly even at very low densities if the temperature is around the 
boiling point. We compare the calculated mass and isotope distributions of 
fragments with the solar abundance [21]. In Sec. 4, we discuss the possibility 
of the ejection of nuclei synthesized in the coexisting region referring to a 
hydrodynamical calculation of supernova explosion [20]. We summarize our 
work in Sec. 5. 
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2 Relativistic Mean Field Approach 



2. 1 RMF Lagrangian and parameter set 

Relativistic Mean Field (RMF) approach has been developed as an effective 
theory to describe the nuclear matter saturation in a simple way [16]. Hav- 
ing improvements to include meson self-coupling terms, it describes well the 
binding energies of neutron-rich unstable nuclei in addition to nuclear matter 
and stable nuclei [17]. 

In this work, we adopt an RMF parameter set TM1 [17]. It has been demon- 
strated that this parameter set TM1 can reproduce nuclear properties includ- 
ing proton- and neutron-rich unstable nuclei. In addition, the EoS table with 
TM1 has been successfully applied to neutron stars and supernova explo- 
sions [13,22-24]. Therefore, it is expected to be reliable also in describing two- 
phase coexistence in supernova matter, which contains asymmetric nuclear 
matter having proton-to-neutron ratio varying in a wide range. 

The Lagrangian contains three meson fields; scalar-isoscalar a meson, vector- 
isoscalar uj meson, and vector-isovector p meson. In this work, we limit the con- 
stituent particles as nucleons, electrons, electron-neutrinos, their anti-particles 
and photons. The explicit form of the Lagrangian including leptons is given 
as follows. 

C =$ N (i0-M- g a a -g u </>- g p r a / a ) ^ N 
+ ^ad^a - 2 m ^ 2 - - -g^a A 

(i0 - m e ) ^ e +$JW V - , 
Wp,, = d p u u - dyU)^ , 

F IM/ = d li A v -d v A li . (1) 

In this section, photon contributions to pressure, energy and entropies are 
included, but we have dropped the photon-couplings to charged particles (See 
Sec. 3 on this point). 

In a mean field approximation, three meson fields are replaced to their 
expectation values. Self-consistency condition for these values can be derived 
in a standard manner, as dP/dx = where x represents the meson field 
expectation value and P denotes the pressure. 
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Supernova matter is composed of nucleons, electrons, neutrinos, their anti- 
particles and photons, characterized by a fixed lepton fraction (lepton-to- 
baryon ratio) Y L , due to the neutrino trapping at high density. As a result, 
there are three conserved quantities; baryon number B, total charge C, and 
lepton number L. Then the chemical potentials of particles are represented by 
the corresponding three chemical potentials, 

Hi = Bifi B + CiHc + Liii L , (2) 

where //#, pc, Pl are baryon, charge and lepton chemical potentials, and Bi, Ci, Li 
are baryon, charge and lepton numbers of particle species i. The chemical equi- 
librium conditions for given conserved densities p k are expressed as dPj dp k = 
Pk- We solve these conditions in a multi-dimensional Newton's method by 
iteration. 

2.2 Two-phase coexistence treatment in RMF 

We apply the mean field approximation to the liquid and gas phases sep- 
arately. In this treatment, we implicitly assume that two coexisting (liquid, 
gas) phases are uniform and have infinite size. We solve chemical and thermal 
equilibrium conditions between the liquid and gas phases. 

In order to make liquid and gas phases coexist at equilibrium, we must 
apply the Gibbs conditions rather than the Maxwell construction, since the 
number of conserved quantity (= 3) is larger than one. The Gibbs conditions 
are given as follows, 

(I - a) p^ + ap° aS = Pk , p L k lq -=p G k a \ pLia.^pGas^ {?>) 

where k = B,C, L. The quantities labeled by Liq. and Gas are those of liquid 
and gas phases, respectively. The gas volume fraction, a, is a number between 
zero and unity. We solve these conditions by using multi-dimensional Newton's 
method, in which the dimension is five in the asymmetric nuclear matter 
(k = B,C) and seven in the supernova matter (k = B,C,L), where the 
variables are p k 19 ' , p k as and a. See Appendix A for the numerical technique at 
very low densities. 

We show the liquid-gas coexisting region of supernova matter in Fig. 1. See 
Appendix B for the method to determine this boundary. The solid curves show 
the boundary of coexisting region, and two phases coexist below the boundary. 
We find that the critical temperature is very high T c ~14 MeV for Yj, = 
(0.3 — 0.4), which is the ratio expected in actual supernova explosions. These 
critical temperatures are much higher than those in neutrino-less supernova 
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Fig. 1. Boundary of liquid-gas coexisting region in supernova matter (solid curves) 
in comparison with symmetric nuclear matter (Y p = 0.5, dotted line) and neu- 
trino-less supernova matter (or finite temperature neutron star matter, m = 0, 
dashed line). 



matter (dashed line) and comparable to those in the symmetric nuclear matter 
(dotted line). The boiling points of supernova matter remain to be Tb oi i ~ 1 
MeV even at very low densities for all lepton fractions. 



The large value of T^oii is due to symmetrization of nuclear matter by 
leptons. Since Yl is kept constant but Y p is not fixed, the supernova mat- 
ter searches its minimum in free energy by changing Y p . In high density 
supernova matter, the lepton fraction is shared by electrons and neutrinos 
(2/?>Y L < Y e < Y L ), then the net neutrino fraction and the lepton chem- 
ical potential become positive (fi L = fi u > 0). Compared to neutrino- less 
supernova matter (fi L = 0), this positive fi L in supernova matter helps to 
enhance Y e {= Y p ) (fi e = —fj, c + /i L ), and to symmetrize nuclear matter. At 
very low baryon densities, electrons become non-degenerate, i.e. electron and 
anti-electron densities are much higher than the net electron densities because 
of small electron mass, and this also applies to neutrinos. Hence, electron and 
neutrino chemical potentials become small, and the charge chemical potential 
also becomes small. This causes nuclear matter symmetric. 

In order to demonstrate this point, we show the density dependence of the 
asymmetry parameter of supernova matter in Fig. 2. Dotted, solid and dashed 
curves show the asymmetry parameter in uniform (homogeneous), two-phase 
coexisting, and the liquid part of coexisting matter, respectively. It is clear 
that, as the baryon density decreases, asymmetry parameter decreases and 
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Fig. 2. Asymmetry parameter 5 = 1 — 21^,, of supernova matter as a function of 
baryon density at T = 1 MeV. Dotted, solid and dashed curves show the asymmetry 
parameter in uniform (homogeneous), two-phase coexisting, and the liquid part of 
coexisting matter, respectively. 

this tendency is stronger in the coexisting region. In uniform matter, the 
asymmetry decreases rapidly at around p B < 10 -7 fm -3 , where leptons domi- 
nate pressure and energy density, as shown later in Fig. 4. This is consistent 
with the above consideration on the lepton dominance. On the other hand, the 
asymmetry parameter start to decrease at much higher density in coexisting 
matter. This symmetrization is mainly due to the symmetry energy in nuclear 
matter. Since the baryon density in the liquid part of matter is around p , 
nucleons can gain symmetry energy by reducing the asymmetry. As shown 
by dashed lines, the asymmetry decreases quickly in the liquid part in the 
coexisting region. 

We turn our attention to one of the characteristic features of the liquid-gas 
coexisting region in supernova matter, which allows many fragments to be 
formed even at very low densities. First we define a new quantity, gas fraction 
Y g , Y g — (baryons in gas phase)/(total baryons) = apg as /p B as a measure of 
bulk fragment yield. Here, we can adopt the normal nuclear matter density po 
for the baryon densities in the liquid phase. After we substitute these baryon 
densities p and p B as for the Gibbs condition, we can obtain the following 
equation, 



Po - Pb 



(4) 



a ~ 



PQ- Pb 
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Fig. 3. Gas fraction Y g in asymmetric nuclear matter (thin line surface) and in su- 
pernova matter (thick lines) as a function of the baryon density and the asymmetry, 
5 = 1 - 2Y p , at T = 1 MeV. 

For example, all baryons are bound in nuclei at Y g = 0. Smaller Y g means 
larger amount of fragments. In Fig. 3, we show the gas fraction in nuclear 
matter composed of neutrons and protons (np matter) and supernova matter 
as a function of the baryon density and the asymmetry parameter, calculated 
in the RMF model. We fix Y p for calculations of np matter. The behavior of Y g 
in np matter is smooth. In symmetric nuclear matter (Y p = 0.5), both of the 
liquid and gas phases are symmetric, and the density in each phase is constant 
in the coexisting region. Then the gas fraction can be expressed by liquid, gas 
and the given average densities {p L B %q ' — po, p B as , p B ) as 



pf as (Po-p B ) 
p B (pv-p G B as ) 



This is a monotonically decreasing function of the baryon density and very 
small in the density range under consideration. When the asymmetry in- 
creases, the liquid phase loses the symmetry energy and nucleons are emitted 
to the gas phase, while Y g is still a decreasing function of p B . 

The behavior of Y g in supernova matter is very different from that in np 
matter. As shown by the thick lines in Fig. 3, as the baryon density decreases, 
Y g first decreases then grows again. This behavior is determined by the proton 
fraction (proton-to-baryon ratio) Y p . In the medium density region (p B = 
(1CT 7 — 1CT 2 ) fm -3 ), matter becomes symmetric as p B decreases, and the gas 
baryon fraction decreases to the symmetric nuclear matter value in Eq. (5). 
Therefore, baryons favor liquid state than nucleon gas in the coexisting region. 
After reaching symmetric matter, gas fraction increases again according to 
Eq. (5). It is also interesting to note that the proton fraction Y p of supernova 
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SN matter EoS (t=o.i,o.5, 1,2, 5, 10, 15,20 Mev) 
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Fig. 4. Equation of state of supernova matter for lepton fractions Yl = 0.1, 0.2, 0.3 
and 0.4. Temperatures are T =0.1, 0.5, 1, 2, 5, 10, 15 and 20 MeV, from top to 
down (from down to top) for F/A (E/A and P/ p B ). Results for uniform supernova 
matter are shown by dotted curves, and those with liquid-gas coexistence are shown 
by solid curves. For energies per baryon {E/A) and free energies per baryon (F/A), 
nucleon mass is subtracted. The pressures are increased by a factor of 10 from low 
to high temperatures (from bottom to top). 

matter shows clear Yl dependence at p B > 1CT 5 fm -3 , while Y v is almost 
independent on Y L at p B < 1CT 5 fm~ 3 . This density roughly corresponds to 
the neutrino- sphere, inside which neutrinos are trapped in supernova core. 

2.3 Supernova Matter Equation of State 

We show the EoS of supernova matter with (without) two-phase coexistence 
by solid (dotted) lines in Fig. 4. In both of the cases, energy and pressure in- 
crease at high densities above p because of the vector meson contributions. 
At very low densities, lepton contributions become dominant, because lepton 
pressure and energy densities are finite even with p B = at finite temper- 
atures. In uniform matter, the nuclear pressure and energy become close to 
free-gas values Pn/Pb - ► T,E N /A — > 3T/2. Between these two extremes of 
density, we can find the effects of two-phase coexistence. 

In two-phase coexistence, by making liquid phase whose baryon density 
is around p , binding energy is gained. At low temperatures, since nuclear 
matter tends to be symmetric and liquid phase is dominant in a wide range 
of densities, liquid part of energy per baryon approaches toward —16 MeV. 
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SN matter Adiabatic Path (S/B=1,2,..., 10,20,100) 
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Fig. 5. Adiabatic paths in supernova matter calculated in RMF. The thin solid 
(dashed) lines show the adiabatic paths of supernova matter with (without) the 
liquid-gas coexistence. Thick solid lines are the boundary of coexisting region. The 
each adiabat corresponds to entropy per baryon S/B = 1, 2, 10, 20, 100. The 
panels are in case of lepton fraction Yl = 0.1, 0.15, 0.5 from upper-left to 
lower-right, respectively. 

At large lepton fractions (Yl = 0.3 — 0.4), the pressure behaves as in the 
case of Maxwell construction for volume instability; pressures in the coexisting 
region are lower (higher) at low (high) densities than in uniform matter during 
the coexistence. On the other hand, we can find double phase transition be- 
havior in lepton deficient matter (Yl = 0.1 — 0.2); while the pressure behaves 
as in the case of volume instability at higher densities (p B > 10 -4 fm -3 , there 
is another overtaking in the density region of p B = 10~ 8 ~ 10~ 4 fm" 3 . At these 
densities, the matter becomes unstable to the small fluctuations of proton frac- 
tion, and liquid part of the matter becomes rapidly symmetric as the density 
decreases. Thus the overtaking of the pressure in the low-density region may 
be suggesting the phase transition in the isospin degrees of freedom. 

2.4 Adiabatic Paths 

After the core bounce, some part of supernova matter expands almost adia- 
batically, and high entropy (S/B > 10) part of the matter would be ejected to 



10 



outside [20]. (See also Section 4.) Therefore, it is important to examine super- 
nova matter along the adiabatic paths. In Fig. 5, we show the adiabatic paths 
in supernova matter with lepton fraction Yl = 0.10,0.15,.. .0.50. Entropy per 
baryon is taken to be S/B — 1, 2, . . . 10, 20, and 100. At very high entropies 
such as S/B = 100, adiabatic paths are almost independent on Y L because of 
the lepton dominance. At lower entropies, we can clearly see nuclear and coex- 
istence effects. While adiabatic paths in uniform matter (dashed lines) evolve 
very smoothly, those with two-phase coexistence (thin solid lines) bend closer 
to the boundary of coexisting region. This can be seen even at S/B = 20. This 
bending comes from the latent heat and suggests that significant amount of 
nuclei are formed around the boundary. 



3 Fragment Distribution in Supernova Matter 

Observations in the previous section tell us that it would be possible for 
supernova matter to experience the liquid-gas phase transition before it is 
ejected to outside. Therefore, it is interesting to study the composition of 
ejecta experienced the phase coexistence. However, we have assumed that the 
coexisting two phases are infinite. It is also to be noted that nuclear matter 
in the liquid phase is more symmetric than in the gas phase. As a result, 
the Coulomb energy is expected to be large, which we have neglected in the 
previous section, and the infinite matter in the liquid phase will fragment into 
finite nuclei. Although the effects of these nuclear formation on EoS may not 
be very large [12], fragment distribution at freeze-out would serve as the initial 
condition for the following r-process. Therefore, in this section, we evaluate 
the fragment yield in a fragment-based statistical model. 

3.1 Statistical Model of Fragments 

In order to describe distribution of finite nuclei, we utilize a fragment-based 
statistical model. This kind of statistical models have been widely used in 
heavy-ion collision studies [4-7] as well as in astrophysics, as referred to the 
nuclear statistical equilibrium (NSE) [11,25]. In NSE, we solve the statistical 
equilibrium condition among nucleons, fragments and leptons in fragment- 
based grand canonical ensemble. In this work, we have ignored relativistic 
corrections and anti-particle contributions of fragments, and fragments are 
assumed to follow the Boltzmann statistics, while leptons are treated rela- 
tivistically. The ensemble averages can be generated from the grand potential, 

Vl = -PV = -VTj^Pf -PtV -PjV , (6) 
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fx f = Z f (fx p -m N ) + N f (p n -m N ) , (8) 



where p/, M/, 5/ and /x/ denote the density, mass, binding energy and the 
chemical potential of fragment /, and Pi and P 7 are the lepton and photon 
pressures, respectively. The intrinsic partition function, C/(T), has been cal- 
culated by using the level density formula for fragments with Af > 5 [6], 

C/(T)=E^ ) exp(-^; W /T) 

i 

oo 

+ ^3 / ^*e-*VT expQyfcF*) , (9) 



A 

A f o 



a/ 



-y (l - c 2^/ 1/3 ) (MeV^ 1 ) , d = 0.2 (MeV" 1 ) , c 2 = 0.8 , 



where = 2jj^ + 1 is the spin degeneracy of the energy level at excitation 
energy E*^ of the fragment species /. 

In NSE, the nuclear binding energy Bf plays an essential role. Fragment 
yields are sensitive to the binding energy modification due to, for example, the 
medium effects in supernova matter. In studies of heavy-ion collisions, since 
the density and its fluctuation are large, the repulsive interfragment Coulomb 
potentials are taken into account explicitly rather than in the form of mass 
modification. In supernova matter under consideration, attractive electron- 
fragment Coulomb potential effects are more important. Since electron den- 
sity is not negligible and almost constant, we ignore inter-fragment Coulomb 
potentials and the electron effects on intrinsic fragment Coulomb energies are 
incorporated in the form of binding energy modification. We have used the 
Wigner-Seitz approximation in evaluating the Coulomb energy correction [12] 
to the binding energy of nuclei adopted in NSE. We assume that the electrons 
are distributed uniformly in a sphere with radius R e f which is determined to 
cancel the charge of the fragment / at a given electron density p e . 



B f ( Pe )=B f (0)-AV f Cou \p e ) , (10) 

Av c oul _ 3 Zje" /3 1 3 \ _ Ro f _ ( Pe \ 1/3 

AVf — 5l^U^-2%> ' ^ = R~r \z^lA f ) ' (U) 

where -8/(0) is the nuclear binding energy in vacuum, and i? / is the nuclear 
radius. 

It is important to note that the Coulomb energy correction, AV^ oul (p e ), 
contains the term proportional to p^ 3 . Because of this functional form, the 
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correction is meaningfully large even at p B = 10~ 6 po- For example, when the 
electron fraction is Y e = 0.3 , the reduction of the Coulomb energy for heavy 
nuclei amounts to 90 % of the total Coulomb energy at p B = p , and the 
reduction is around 10 MeV even at 10~ 6 po- The increase of binding energy 
acts to enhance heavy nuclei, and some nuclei beyond the dripline at vacuum 
or unstable against fission, can be stabilized in supernova matter. Finite gas 
nucleon density also plays a role to form nuclei beyond the dripline tentatively 
by the balance of nucleon absorption and emission. We have adopted the mass 
table of Myers and Swiatecki [26] , which is based on the Thomas- Fermi model 
with shell correction for about 9000 kinds of nuclei. 

Since nuclear binding energies depend on the electron density, we have to 
solve the chemical equilibrium condition of nuclei and leptons in supernova 
matter in a consistent way to satisfy = p p + p e — p n — p v = 0. Provided 
that the baryon density and temperature are given, and that the charge and 
lepton densities are fixed as (p c ,p L ) = (0, Y^ps), all the particle densities are 
determined if the average proton fraction (Y p ) is given, 

(1 - Y p ) p B = J2 N fPf(Pn, Pp, B f (p e )) = p n (Pn, Pp, Pe) , (12a) 

/ 

Y vPB = Y, Z fPAP^PpiB f (p e )) = p p (p n ,p p ,p e ) , (12b) 
/ 

Y p p B = p e (p e ) , (Y L - Y p ) p B = pu(Pu) , (12c) 

where p n and p p are the neutron and proton densities including those in nu- 
clei. We can easily solve last two conditions in Eqs. (12c) numerically, and 
derivatives of p e and p u with respect to Y p can be obtained as dp e /dY p = 
p B (dpe/dpe) -1 and dp u /dY p = —p B idp v j 'dp u )~ l . Once Y p is given, first two 
equations are the same as usual conditions in fragment-based statistical mod- 
els. By using the charge neutrality condition p e = Y p p B , we can get the deriva- 
tives of p„ and p p with respect to Y p as, 



^ dpjdpn , dpjdp p ^ 




-l-d-pjdp e ^ 



= p B dY p 

\ dpp/dpn , dpp/dpp J \dp p J y 1 -dp p /dp e J 



(13) 



Therefore, we can solve the chemical equilibrium condition, F M = 0, by the 
Newton's method, 5F^ = -Fj {dF^ (Y p ) jdY p ). 



3.2 Gas Fraction and Equation of State in a Statistical Model 



Compared to RMF, there is no sharp phase transition in NSE because of the 
finite size of fragments. Although we can see kink-like behavior in the nuclear 
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Fig. 6. Left: Mass moment ratio (upper panel) and caloric curve (lower) of supernova 
matter with Yl = 0.35 at p B = 10~ 7 and 10~ 3 fm~ 3 . In the upper panel, mass 
moment ratios, M 2 /Mi (thin solid lines), £ = a 2 (Af)/ < Aj > (with a, dotted 
lines), £' (same as £ but without a, thick solid lines) are shown. Right: £' (upper 
panel) and caloric curve (lower) of supernova matter with Yl = 0.35 at baryon 
densities p B = 10~ 7 , 10~ 5 , 10 -3 and 10 -2 fm -3 . Results with Coulomb correction 
(solid lines) and without Coulomb correction (dashed lines) are compared. 
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Fig. 7. Liquid-gas coexisting region of supernova matter with Yl = 0.35 calculated 
in NSE (solid line) in comparison with RMF results (dashed line). Results without 
Coulomb correction is shown by dotted curve. Boiling point in the statistical model 
has been determined as the maximum point of mass variance-to-average ratio, while 
there appear two local maxima in the fragment mass fluctuation when we include 
the contribution of a particle, as shown by dot-dashed curves. 

part of energy per baryon as a function of the temperature (caloric curve), 
this kink is not clear enough to define the boiling point as seen in Fig. 6 . 
There are several definitions of T hoi i proposed in the literature. For example, 
T boi i is proposed to be well defined at the peak of M 2 /M 1 by Bauer [7]. The 
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Table 1 

Boiling points in NSE at Yl = 0.35 as a function of the baryon density. Results 
with (NSE) and without (NSE, nc) Coulomb corrections are shown. 



p s (fm- 3 ) 


10- 7 


10" 5 


10~ 3 


io- 2 


T 6oji (NSE) (MeV) 


0.77 


1.43 


3.90 


6.88 


T ftoii (NSE,nc) (MeV) 


0.77 


1.40 


3.52 


5.42 



n-th moment of light fragments, M n , is denned as M n = S/^/P/- Another 
way is to define T^u by the peak of the light fragment mass variance-to- 
average ratio, £ = a 2 (Af)/ < Af >, which becomes small when one fragment 
(or nucleon) dominates, and becomes unity when the mass distribution is a 
Poisonnian. We show these ratios in the upper-left panel of Fig. 6. At high 
densities, these definitions give reasonable boiling points, but at low densities, 
the peak of M 2 /M 1 becomes dull, and £ shows two peaks. From fragment 
distributions, we find that the two peak structure appears due to the formation 
of a which can be comparable to neutrons at very low densities, where gas 
part becomes almost symmetric. In Fig. 7, we show the temperatures at local 
maxima of £ by dot-dashed lines. This formation of a makes also MijM\ 
peak dull. Therefore, we here define the boiling point as the peak of the light 
fragment mass variance-to-average ratio where we excluded a particle in 
the calculation of the light fragment mass average and variance. As shown in 
the solid lines in Fig. 6, the peak position of £' is well defined at any density. 
In addition, this boiling point corresponds to the kink position in the caloric 
curve. 

In Fig. 7, we show the density dependence of T^n, in comparison with the 
RMF results. We find that the boiling points in NSE are lower than those in the 
two phase treatment of the RMF model by about a factor of two. The reduction 
of Thou is a natural consequence of finite size of nuclei. Because the intrinsic 
Coulomb energy cannot be completely removed by electrons at densities p B < 
Po, nuclei are limited to have finite size. Then nuclei loses surface energy in 
addition to the Coulomb energy, and gas nucleons are favored. However, it is 
worthwhile to note that the qualitative behavior of T^u is similar to that in 
RMF, and they are still high enough, T boU > 0.7 MeV for p B > 10~ 7 far 3 . We 
tabulate the boiling points at p B = 10~ 7 , 10~ 5 , 10~ 3 and 10~ 2 fm -3 in Table 1. 

fa the RMF treatment, one of the most characteristic features in the coex- 
isting region is the reduction of the gas fraction. This also applies to the NSE 
results. We have defined the gas fraction as the ratio of isolated nucleon density 
to the total baryon density, Y g = (p p + p n )/ p B . As shown in Fig. 8, gas fraction 
behaves similarly to that in RMF results; as the baryon density decreases, it 
quickly becomes very small until p B ~ 10~ 8 fm~ 3 , and gradually grows at lower 
densities again, fa addition, it is interesting to note that Y g curves within the 
statistical model have the second minimum at p B ~ 10~ 3 fm~ 3 . Supernova 
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matter favors nuclear fragment (nucleus) state rather than nucleon gas at 
these densities. We call these minimum regions as the first (p B ~ 10~ 7 fm -3 ) 
and the second (p B ~ 1CT 3 fm~ 3 ) fragment windows, respectively. The first 
one is caused by the drastic isospin-symmetrization of supernova matter by 
leptons at low baryon densities. The second fragment window is specific to 
NSE. The mechanism of this appearance is not very clear, but we find that 
these two minima converges to one when we use the liquid drop mass formula 
for the nuclear binding energies and ignore the surface term. 

Y L =0.35, T=0.5,0.6,.. .,1.0,2.0 MeV 

< n 
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Fig. 8. Gas fraction in supernova matter as a function of the baryon density at 
Yl = 0.35. Results in NSE (RMF) are shown by solid (dashed) lines. Temperatures 
are taken to be T = 0.5, 0.6, 1.0 and 2.0 MeV from bottom to top. 

In Fig. 9, we show the EoS in NSE at T = 5 and 1 MeV with Y L = 0.35, in 
comparison with those in RMF with the Thomas-Fermi approximation for a 
dominant configuration [13] (RMF+TF, dotted lines), RMF with two-phase 
coexistence (RMF(coex.), dashed lines), and homogeneous RMF (RMF(unif.), 
dot-dashed lines). We find that the finiteness of nuclei does not affect the EoS 
at high densities, p B > 10~ 2 fm" 3 (10~ 4 fm" 3 ) for T = 5 MeV (1 MeV), but 
modifies the density dependence of the pressure at lower densities, as seen in 
the difference between RMF(coex.) and RMF+TF. Compared with the EoS 
in RMF+TF, the present NSE results give very similar pressures, except for 
the density region p B ~ 10~ 3 fm+ 3 (10~ 6 fm~ 3 ) for T = 5 MeV (1 MeV), where 
the pressure are different by 10 ~ 20 %. These densities correspond to the 
region where the given temperatures are close to the boiling points. 

The agreement of EoS in NSE and RMF+TF is somewhat surprising. There 
are three large differences in NSE and RMF+TF: (1) Nuclear masses are taken 
from the table [26] in NSE, while masses are calculated in RMF+TF. (2) In- 
terfragment and nucleon-fragment nuclear interactions are neglected in NSE, 
while they are included in RMF+TF. (3) One configuration is assumed in 
RMF+TF, while statistical ensemble is considered in NSE. Thus the above 
agreement might suggest that once nuclear masses are properly included, nu- 
clear interactions between gas nucleon and fragments play a minor role in 
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Fig. 9. Equation of state of the supernova matter in NSE in comparison with those 
in RMF models at T = 1 MeV (lower curves) and T = 5 MeV (upper curves, 
scaled up by a factor of 10). The solid, dotted, dashed and dot-dashed line denote 
EoS in NSE, RMF with Thomas-Fermi approximation model (RMF+TF) [13], the 
liquid-gas coexisting RMF model (RMF(Coex.)) and homogeneous RMF model 
(RMF(Unif.)), respectively. 



3.3 Fragment Distribution and Coulomb Correction Effects 



In this subsection, we investigate fragment distribution in the density range 
10~ 7 < p B < 10~ 2 fm -3 , starting from the first fragment window to the density 
close to the critical point. We are most interested in the boundary of coexisting 
region at the temperatures around Tb oi i(p B ). If the density is not very small, 
the freeze-out temperature Tf is expected to be around or just below the 
boiling point T^ou- Below the boiling points, fragments are formed abundantly 
by absorbing many of gas nucleons, and nuclear number density of particles 
(sum of nucleon and fragment densities) becomes quickly smaller, then the 
mean free path for each particle becomes much longer. This rapid fragment 
formation makes the typical interaction intervals longer, and is expected to 
help the system to freeze-out. The condition of freeze-out should be studied 
more carefully. 

We can see common features of the temperature dependence of fragment 
mass distributions in each row of Fig. 10. When the temperature is higher than 
the boiling point, the fragments almost obey the exponential distribution. At 
the temperature near the boiling point, the distribution shows the power-law 
like behavior up to some mass. A similar trend, Y{Af) ~ ^4/ T ) a l so appears 
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in nuclear multifragmentation [1], where Y and Af denote the fragment yield 
and mass. This power law was suggested by Fisher [2] for a mass distribution 
of droplets at around the critical point. As the temperature becomes lower 
than the boiling point, the distribution at each density becomes localized to 
some mass number. The localization is caused by the small entropy and shell 
effects. 



p B = 10" 7 fm" 3 p B =10" 5 fm" 3 p B = 10" 3 fm" 3 p B = 10" 2 fm 

(1st Frag. Win.) (v-sphere surface) (2nd Frag. Win.) 
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Fig. 10. The fragment distribution at p B = 10~ 7 fm~ 3 (first column), 1CP 5 fm -3 
(second column), 10~ 3 fm~ 3 (third column), and 10~ 2 fm~ 3 (fourth column) with 
lepton fraction Yj, = 0.35 in comparison with the observed solar abundance (dots). 
Temperatures are taken to be T = 0.2 Tt> u ~ 1.1 Tf, n, where Tf, n is the boiling 
point at each density. Thick and thin lines denote the distribution with and without 
Coulomb correction, respectively. 

The density p B = 1CT 7 fm~ 3 corresponds to the the first fragment window, 
where supernova matter is symmetrized by leptons and nuclei are formed 



18 



abundantly at low temperatures. As shown in Fig. 10 (first column), the main 
products near the boiling point (T = T^i ~ 0.77 MeV) are nucleon, a, iron 
peak nuclei, and they extend to the first peak nuclei of r-process in the dis- 
tribution shoulder. The Coulomb correction effect is negligible at this density. 
Although the above distribution seems to be similar to the initial seed nuclear 
composition in a standard scenario of the r-process, we would like to point out 
that when the temperature rises to around T^n ~ 0.8 MeV, the distribution 
become broad due to thermal fluctuation. Especially, just above the boiling 
point, T = 1.1 T boi i (top panel), the equilibrium nuclear distribution resembles 
to the solar abundance up to the iron peak. 

Next, we show the results at p B = 10~ 5 fm~ 3 in Fig. 10 (second column), 
which corresponds to the density around the neutrino sphere. Most stable 
nuclei at this density are those at around the first peak of r-process, as can 
be seen in the distribution at low temperatures. As a result, these nuclei are 
formed easily also at around the boiling point (T = 1.0 Tj, n = 1.43 MeV), 
in addition to nucleons, a, and iron peak nuclei. The Coulomb correction 
effect is not large, but at around T hoih we can see small enhancement of heavy 
fragments with larger masses over the first peak of r-process. 

In the second fragment window (p B = 10~ 3 fm~ 3 ), power-law like behavior 
can be seen at around Tj, u = 3.90 MeV up to the second peak of r-process, as 
shown in the third column of Fig. 10. Coulomb correction effects are clearly 
seen at this density. The center of the distribution shifts from the first peak 
of r-process without Coulomb correction to the second peak of r-process with 
Coulomb correction. 

As shown in the fourth column of Fig. 10, the distribution at p B = 10~ 2 fm -3 
shows the same trend. It is interesting to find that the peak position in the NSE 
results is shifted downwards compared to the observed third peak of r-process. 
This shift mainly comes from the np ratio of formed nuclei. The observed third 
peak of r-process is a consequence of the neutron magic number N = 126 and 
the np ratio along the r-process path. In the present NSE model calculation, 
having a large np ratio, nuclei beyond the neutron dripline appear easily at 
equilibrium. 

We also give an example of isotope distribution at Yl = 0.35 in Fig. 11. Here, 
we choose a slightly lower temperature than the boiling point, Tf = 0.9 T^, 
for the freeze-out temperature, as discussed at the beginning of this subsection. 
The upper- left, upper-right, lower- left and lower-right panels of Fig. 11 show 
the isotope distributions at p B = 10~ 2 , 10~ 3 , 10~ 5 and 10~ 7 fm~ 3 , respectively. 
At p B = 10~ 2 fm~ 3 , while the observed isotope ratio are well explained in the 
calculation from Z = 20 (Ca) to Z = 92 (U) with one overall normalization 
factor, much more neutron rich nuclei appear at equilibrium. For isotones with 
N = 126, nuclei with Z = 54 to Z = 92 are formed. This large np ratio, which 
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comes from the large np asymmetry of the liquid phase as shown in Fig. 2, 
gives smaller mass number with N = 126. 



At lower densities, similar trends can be seen in the mass number range 
produced at around Tboii, except for p B = 1CT 7 fm~ 3 . As already discussed, 
nuclear matter is symmetrized outside of the neutrino sphere, then the calcu- 
lated distribution at p B = 10~ 7 fm~ 3 is shifted toward proton rich side of the 
observed solar abundance for large Z. 
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Fig. 11. Calculated even Z isotope distribution at p B = 10~ 2 , 10 -3 , 10 -5 and 
10~ 7 fm~ 3 and (T, Yl) = (0.9 Tj, oi i, 0.35) in comparison with the solar abundance by 
circle and triangle symbols. The boiling points are those for given densities. Yields 
are shifted, for clarify of plot, by a factor of 10 from Z = 20 (Ca) to larger Z nuclei 
(from top to bottom). One overall factor is chosen to get good global fit. 



3.4 Fragment Distribution outside of Neutrino Sphere and in Neutrino-less 
Supernova Matter 



While we have considered the /9-equilibrium condition with fixed lepton frac- 
tion Y L , proton fraction Y p is almost independent on Y L outside of the neutrino 
sphere. Once Y p , p B , and T are given, neutron and proton chemical potentials 
are uniquely determined, thus boiling point and fragment distribution are also 
obtained uniquely. In RMF, this independence on Y L outside of the neutrino 
sphere can be found in Y p (Fig. 2) and in T boi i (Fig. 1). The same tendency 
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Mass Dist. in SN Matter,(p B , T)=(10" 7 fm" 3 , 0.9T boN ) 
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Fig. 12. The fragment distribution in neutrino-less supernova matter with 
Y p = 0.46 (thick lines, n,p,e~) and supernova matter (thin lines, n,p, e~,v e ) at 
p B = 10~ 7 fm -3 in comparison with the solar abundance (filled circles). For su- 
pernova matter, we have chosen Yl = 0.1 ~ 0.45. The boiling points are almost 
independent on Y L , T boU ~ 0.77 MeV. 

applies to fragment distribution in the NSE results of supernova matter and 
neutrino-less supernova (NS) matter in which neutrino chemical potential is 
zero (fi u — /i L — 0). In Fig. 12, we show the fragment mass distribution in 
NS matter at p B = 10~ 7 fm" 3 , in which Y p = 0.46, in comparison with that 
in supernova matter (Y L = 0.1 ~ 0.45). We find that the boiling points and 
nuclear distribution are almost independent on Y L , and they are almost the 
same as those in NS matter. 

The above observations tell us that the fragment distributions under ther- 
mal and charge equilibrium at around T hoa are independent of degree of neu- 
trino trapping. In addition, the consequent fragment distribution provides 
iron-group nuclei, which become seed elements in a standard r-process sce- 
nario. On the contrary, Y p clearly depends on Yl inside the neutrino sphere. 
The nuclear distribution formed in the neutrino sphere is sensitive to the dy- 
namics of supernova explosion, especially on how much neutrinos escape before 
the neutrino trapping. 



4 Relation to Nucleosynthesis in Supernovae 

One of the key questions on the LG-process is whether the elements made 
through the liquid-gas phase transition are ejected to outside. If they are 
ejected, they contribute to the following nucleosynthesis such as the r-process. 

As mentioned earlier, all trajectories of adiabatic path go across the bound- 
ary of coexisting region. If the material with S/B > 10 is ejected in super- 
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nova explosion, the present model calculation suggests that the material would 
go through the coexisting region before the ejection. This can happen if the 
freeze-out temperature of supernova matter is as low as T = (0.5 — 2) MeV 
at densities p B = (1CT 7 — 1CT 5 ) fm~ 3 , where S/B > 8 at the boiling points 
provided that the adiabatic paths go straight until the boundary. Since the 
fragment distribution is almost independent of the lepton fraction at these 
densities, the distribution shown in Fig. 10 would appear regardless of details 
of neutrino trapping. 

In order to examine this possibility, we analyse the results of hydrodynam- 
ical calculation of core-collapse supernova [20]. In their calculation, the EoS 
table derived by the RMF-TF model with the TM1 interaction [13] has been 
used in the general relativistic hydrodynamics [27]. Note that this EoS table 
and the EoS in NSE give almost the same pressure as shown in Fig. 9. Among 
series of model calculations of adiabatic collapse, a case of iron core of 15M 
presupernova star [28] explode hydrodynamically with the fixed electron frac- 
tion, which were assumed to study the properties of EoS in model explosion. 
In this model explosion, the material from deep inside the core can be ejected. 
The entropy per baryon of this ejecta is about S/B ~ 10 determined by the 
shock passage. The trajectory of this ejecta passes through the coexisting re- 
gion. Further inner material has lower entropies down to S/B ~ 1, which 
are favorable for the LG process, but the ejection becomes rather difficult. It 
is also to be noted that the model explosion is obtained by a simplified cal- 
culation without neutrino-transfer assuming large electron fraction. Further 
careful studies are necessary to examine quantitatively the mass ejection by 
performing the numerical simulations of hydrodynamics with neutrino-transfer 
of core-collapse. 

Although hydrodynamical explosion (so-called prompt explosion) has been 
claimed to be limited for very small iron core and extreme cases of EoS and 
electron capture rates [8,9], the outcome of nucleosynthesis is extremely in- 
teresting if the mass ejection really occurs. A case of prompt explosion of 
the small iron core of 11M© presupernova star has been studied as a site 
of successful r-process nucleosynthesis by the dynamical ejection of neutron- 
rich material [29]. In their studies, the material having entropy per baryon 
S/B ~ 10 is ejected. These ejecta may be affected by the LG-process to help 
the production of heavier r-process elements around actinides. If the ejection 
from relatively high density (p B > 10~ 5 fm~ 3 ) takes place, the formed frag- 
ments are generally very neutron rich and some part of them are unstable 
against neutron emission. These nuclei provide a huge amount of neutrons, 
which help the following r-process to proceed. Mass ejection is also expected 
in asymmetric supernova explosion. The convection of material is generally 
seen in supernova simulations and the material deep inside may be ejected by 
hydrodynamical instability. (See [30] for example.) Jet-like ejection of material 
may occur in rotating collapse of stars as discussed in some literatures [31,32]. 
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The ejection of material and its consequence on nucleosynthesis should be 
clarified together with the problem of supernova explosion mechanism and 
require further extensive studies. 



5 Summary and discussion 

In this paper, we have investigated the liquid-gas phase transition of su- 
pernova matter, and its effects on the fragment formation. We have used two 
models — the Relativistic Mean Field (RMF) model and the Nuclear Statis- 
tical Equilibrium (NSE) model. 

In RMF, we have used the interaction TM1, which has been successfully 
applied to finite nuclei including neutron rich unstable nuclei, neutron stars, 
and supernova explosion [13,22-24]. Leptons are shown to play non-trivial 
roles such as the symmetrization of nuclear part of supernova matter. As a 
result, nuclear liquid gains symmetry energy, and the calculated boiling points 
in supernova matter {T^i > 1 MeV for p B > 10~ 10 fm~ 3 ) are comparable 
to those in symmetric nuclear matter at low densities. Adiabatic paths are 
shown to go across the boundary of coexisting region even at high entropy 
such as S/B > 10, which is expected to be enough for supernova matter to be 
ejected to outside. Clear concentration of adiabatic paths to the boundary of 
coexisting region have been found. All of these findings suggest that at least 
a part of ejecta in supernova explosion would experience the liquid-gas phase 
transition before freeze-out. 

In NSE, we have used nuclear binding energies of Myers- Swiatecki model [26] 
with Coulomb correction due to electron screening as a medium effect [12]. 
Since larger species of nuclei become stable with this Coulomb energy correc- 
tion, we have adopted the mass table of around 9000 nuclei constructed by 
Myers and Swiatecki [26]. Because of the finiteness of nuclei, they lose surface 
and Coulomb energy compared to the case of coexistence treatment of two infi- 
nite matter phases in RMF. The boiling points become slightly lower, but they 
are still high; T boi i > 0.7 MeV for p B > 10~ 7 fm~ 3 . Calculated fragment mass 
distributions around Tt, oi i(p B ) show enhancement of the iron peak elements, the 
first, second, and third peak r-process elements at p B = 10~ 7 , 10~ 5 , 10~ 3 and 
10~ 2 fm~ 3 , respectively. In addition, calculated isotope distribution shows that 
very neutron rich nuclei around and beyond the neutron dripline may exist 
under thermal and chemical equilibrium in supernova matter with degenerate 
neutrinos. These unstable nuclei against neutron emission would provide a lot 
of neutrons after freeze-out, which may help the r-process to proceed. 

From the present investigations, we can draw a new scenario for making 
seed nuclei before the r-process; fragments are abundantly formed through 
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the liquid-gas phase transition of supernova matter before the freeze-out, and 
this formation of fragments serve to produce the bulk structure of the seed ele- 
ments. We call this process as the LG process as a pre-process of r-process [15]. 

It is interesting to note that our model based on the liquid-gas coexisting 
state of supernova matter can even provide the r-process nuclei or their seed in 
a simple manner based on the condition determined by the dynamics of super- 
nova explosion such as p B , Yl, and T. One of the most promising conditions 
is p B = 1CT 5 fm~ 3 . This density roughly corresponds to the neutrino sphere. 
The entropy at is a little smaller than the ejection criteria, S/B > 10 
in one-dimensional hydrodynamical calculation of supernova explosion [20]. 
However, it would be possible that matter with small entropy can be ejected 
by convection and/or jet in asymmetric supernova explosion [30,32]. The most 
conservative freeze-out density for ejection would be p B = 10~ 7 fm -3 . The en- 
tropy at T boi i is large enough, and the seed nuclei will be nucleons, a, iron 
peak nuclei and a small amount of the first peak nuclei of r-process. Higher 
densities may not be relevant to ejection, but it may be closely related to the 
nuclear distribution on hot neutron star surface. 

In this work, we have assumed equilibrium throughout this paper. One of 
the key questions is the freeze-out conditions of supernova matter, at which 
nuclear reactions become less frequent and supernova matter goes off equilib- 
rium in the expansion time-scale. The seed nuclear distribution of the r-process 
will be given as the nuclear distribution on the freeze-out line in the (p B ,T) 
diagram. It is important to determine the freeze-out condition in supernova 
dynamics. Another important direction is to construct a model which includes 
both of the mean field nature such as in RMF and the statistical nature in 
NSE. In a present NSE treatment, only the Coulomb correction is included as 
the medium effects, and medium effects from strong interactions are neglected. 
This neglection may lead to the overestimate of neutron rich nuclei, as dis- 
cussed in recent statistical fragmentation models [33]. On the other hand, in 
the Thomas-Fermi treatment of heavy-nuclei with EoS derived using RMF, 
since statistical nature or fragment distribution is not taken care of, the treat- 
ment is not sufficient especially at around T hoa . Works in these directions are 
in progress. 
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A Low Density Approximation 

In the calculation of liquid-gas coexistence at low temperatures, it becomes 
necessary to solve the chemical equilibrium of liquid and very low density gas. 
For example, the gas baryon density becomes around p B as = 10 -74 fm -3 at 
the coexisting condition (p s , T) = (10~ 10 fm~ 3 , O.lMeV). In order to efficiently 
obtain the derivative matrix dpi/dpj = B,C,L), which are required in 
solving coexistence, we take the low temperature and low density approxima- 
tion for nucleons at low baryon densities (p B <C 1CT 5 fm~ 3 ) around and below 
the boiling point in the mean field calculation. 

In the mean field approximation, nucleon distribution is a function of the 
effecitve mass M* = M + g a a and effective chemical potentials z/j = pi — 
Quj^q ~ 9 P T iP3o (i = n or p), where o" ,cj and p 30 are the expectation values 
of the meson fields of a, uo and the neutral p mesons, respectively. For a given 
baryon density p B and proton fraction Y p , the vector meson expectation values 
are uniquely deterimed as ujq = ou (p B ),p 30 = p3o(p T ), where p T = p p — p n . In 
a non-relativistic limit [18], we can take energy as E* = M* +p 2 /2M*. At low 
densities (p B <C 10~ 5 fm~ 3 ), the nucleon fugacity fa = exp ( — (M* — v^/T) 
becomes much smaller than unity, then we can safely ignore the second and 
higher order terms in the fugacity f{. In this approximation, the Fermi dis- 
tribution is approximated to be the Boltzmann distribution, then the baryon 
number density pi and the scalar density p s are analytically obtained for a 
given value of a as, 



p s K,^o) =G + {y n , M> ),T) + G+{y p , M*(a ),T) + 0{f) , (A.2) 



The self-consistent condition a = cr (p s (i> n , v p , <r )) can be solved by interation 
using (A.1-A.3), which converges in a few steps at low densities. All of the 
above densities are represented by three variables, v n , v p and a for a given 
temperature T, and we can eliminate the cr dependence by using the total 
derivative of the above self-consistent condition. 



Pi {yi, <7 ) = G.(vi, M*(a ),T) + 0(f?) 
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By solving the first equation (A. 4) in dao and substituting it in the second 
equation (A. 5), we obtain dpi/dfj and then dvi/dpj. Finally, the third equa- 
tion (A. 6) gives the partial derivatives dpi/dpj by eliminating du. Having 
the (relativistic) lepton integrals and the above derivatives in nucleons, it is 
straightforward to construct the matrix dpi/dpj = B,C,L). 



B Procedures to obtain coexisting region of supernova matter in 
RMF 

The boundary of coexisting region of supernova matter has been determined 
in three steps; symmetric nuclear matter, asymmetric nuclear matter, and 
supernova matter. 

First, for a given temperature, we solve the coexisting condition in sym- 
metric nuclear matter {Y p tq ' = Y p as = 0.5). If there is a density region where 
pressure is decreasing for increasing p B , liquid and gas phases can coexist, 
and we can find coexisting densities, p B /iq ' and Pg as , by using the Maxwell 
construction. Secondly, we solve the coexisting condition for (p B iq ', Y p Liq -) and 
(p^ as , Y p Gas ) in asymmetric nuclear matter. Having the coexisting condition at 
a given Y p n \ it is easy to find coexisting condition for slightly different Y p %q ' 
by using the multi-dimensional Newton's method. Starting from symmetric 
nuclear matter, three variables (pg iq ', Pg as , Y p Gas ) are determined for a given 
Y p Ltq - which is slightly different from that in the previously solved condition 
by requiring the condition, pj^ q - = p Gas , pj£ q ' = p Gas , and P Ltq ' = P Gas . 
We show the boundary of coexisting region of (p B , Y p ) by the thick solid line 
in Fig. B.l. Filled circles show the point where two phase become uniform, 
(Pg tq ', Yp tq -) = (p Gas , Y Gas ). The density gap p^ iq - — p Gas generally decreases as 
the liquid becomes more asymmetric, because of the symmetry energy loss. At 
lower temperature, the coexisting pressure and thus the coexisting gas baryon 
density become small, then larger density gap appears. Thirdly, proton frac- 
tion Y p is determined as a function of p B in uniform (homogeneous) supernova 
matter at a given Y L , by using the charge neutrality condition, p e = p p = Y p p B , 
and the chemical equilibrium condition, p u = p L = p e + p p — p n . When Y p is 
in the coexisting region of nuclear matter, liquid and gas phases can coexist 
in supernova matter. The boundary of coexisting region for a given T is deter- 
mined by the crossing point of these two lines. Since Y p becomes smaller for 
smaller Yl as shown in Fig. B.l, the coexisting density region becomes nar- 
rower for smaller Y L . This is the reason why the boiling temperatures decrease 
for smaller lepton fraction as shown in Fig. 1. 

For neutrino-less supernova matter (NS), the procedure is almost the same, 
except that the chemical equilibrium condition is modified to p v = 0. 
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Fig. B.l. The boundary of liquid-gas coexisting region in nuclear matter (thick 
solid lines), and proton fraction as a function of the baryon density in supernova 
matter without (thin solid lines) and with (dotted line) coexistence. Filled circles 
show the points where liquid and gas phases converges to the uniform matter. 
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